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spacing). We propose methods to separately calculate these LECs, typically called Wg, and 
Wg. While only one linear combination of these constants enters into physical quantities, different 
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pion scattering, which are accessible only in the partially quenched extension of the theory. 
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I. INTRODUCTION AND SUMMARY 

Calculations using (improved versions) of Wilson 
fermions [T] have successfully approached P], and even 
reached |3], physical light quark masses. The twisted- 
mass extension has also been highly successful [1] . Never- 
theless, the explicit breaking of chiral symmetry can lead 
to significant lattice artefacts that need to be understood 
and controlled. For example, studying the long-distance 
behavior of Wilson fermions using chiral effective the- 
ory [51 one finds that, when quark masses satisfy 
771 ^ o^Aqqy) [o- being the lattice spacing), discretiza- 
tion errors lead to a non-trivial phase diagram, with one 
scenario (the "first-order scenario") having a minimum 
pion mass, M™™ ~ qAq^y^j ^-nd the other having a re- 
gion of Aoki phase in which flavor is spontaneously bro- 
ken [7j. There have also been numerous studies of the 
properties of mesons and baryons using the chiral ef- 
fective theory — usually called Wilson chiral perturbation 
theory (WChPT) — which provide the functional forms 
needed to do simultaneous chiral and continuum extrap- 
olations Q 

In the unquenched theory with two light flavors (a 
class of theories which includes physical QCD if one 
treats the strange quark as heavy) , the chiral Lagrangian 
contains only one independent term proportional to a^. 
As a result, corrections enter with a single low en- 
ergy constant (LEC), denoted C2 in Ref. 5 . The sign 
of this constant determines the vacuum structure when 
m ~ o^Aq^^. If one wants to go beyond the phase struc- 
ture, however, and use the chiral effective theory to de- 
termine discretization errors in the spectral density of the 
Hermitian Wilson-Dirac operator ^ , or, more generally, 
to determine the detailed properties of low lying eigen- 
values of the Wilson-Dirac operator (THl E] , then it turns 
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FIG. 1. Quark contractions contributing to n n scattering. 
Contractions involving the interchange of the final state pious 
are not shown. 



out that one must consider the partially quenched (PQ) 
extension of WChPT. In this extension, there are three 
LECs entering in terms proportional to a^, denoted Wq, 
W'j and W'^ [see Eq. ( [T5| below], of which ci is a par- 
ticular linear combination [see Eq. (16l]. The detailed 
properties of the spectrum and eigenvalues depend on all 
three LECs and not just on c^- It is thus of interest to 
determine the LECs separately, and this is the topic of 
the present work. 

Our proposal builds upon one of the methods be- 
ing used to determine C2. This is to calculate cer- 
tain pion scattering lengths (i.e. the scattering ampli- 
tudes at threshold), which in the continuum are propor- 
tional to m, but which also have contributions propor- 
tional to c^a^ when discretization errors are included |12j . 
Computationally, the simplest choice is to calculate the 
7r+7r+ —5- 7r+7r+ scattering amplitude, for this involves no 
quark-antiquark annihilation contractions, as illustrated 
in Fig. [T] The scattering length can be calculated from 
the energy shift 5E = i?(7r+7r+) — 2r7i^ using the method 
of Ref. [13] ■ (More details of this method will be given 
below.) Such a calculation is in fact presently being car- 
ried out [n]. 

Our method for calculating two of the LECs is to 
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FIG. 2. Ratios of partially quenched correlation functions 
used to extract LECs. Lines show quark propagators, with 
the valence flavor indicated by the label [and color]. All in- 
terpolating operators are at zero three-momentum, and are 
placed at the Euclidean times indicated. Expectation values 
are taken with two sea quarks having the same action and 
masses as the valence quarks. 



that contributions from excited pion states, which fall 
exponentially, can be neglected. It also requires that 
the linear terms in t can be distinguished from quadratic 
terms which appear at higher order in ChPT and which 
scale as In practice, these conditions seem reason- 

able (see, e.g., Ref. jMj). 

In a ratio of physical correlation functions, the contri- 
bution linear in t is the first term in the expansion of 
the exponential exp{—6Et). In the PQ theory, by con- 
trast, neither nor Rs{t) are exponentials. This 
is immediately clear for Rs{t) due to the absence of an 
i-independent constant term and will be shown explic- 
itly for Rvit) by calculating the quadratic term — see 
Eq. (104). Nevertheless, as long as it is possible to pick 



out the term linear in t one can extract the LECs. 



consider separately the quark-disconnected contraction 
of Fig. IlTa) and the quark-connected contraction of 
Fig. [lJb)T^ This separation is simple in a numerical sim- 
ulation, but it introduces a problem in the theoretical 
description. To pick out the separate contractions re- 
quires (for the two-flavor theory that we consider) using 
a PQ theory. We stress that in this application of partial 
quenching the valence quarks are all degenerate with the 
sea quarks and have the same action. This is in contrast 
with the most common application in which the valence 
and sea quark masses (and sometimes also actions) differ. 

While well defined as a Euclidean field theory, the PQ 
theory is unphysical, and the method of Ref. [13 for de- 
termining the scattering lengths does not apply. In par- 
ticular, the individual contractions cannot be written as 
a sum of exponentials, as is the case for their sum. Our 
proposal is instead to directly fit the correlation func- 
tions calculated in the simulation to the predictions of 
PQWChPT. These predictions depend, at leading order 
(LO), on the LECs Wq and W^, and thus these two con- 
stants can be determined from lattice data at sufficiently 
small m and a^. This method of comparison was first 
introduced in the quenched theory [15] . In both the 
quenched and PQ cases, the key point is that the ef- 
fective chiral theory reproduces, at long distances, the 
unphysical nature of the underlying theory. 

Specifically, we suggest calculating the ratios of corre- 
lation functions shown in Fig. [2] [and defined in Eqs. (51 ■ 
551 below]. This calculation must be done outside the 



Xoki phase, if present, as we assume that flavor is not 
spontaneously broken. As we show in Sec. |IV[ at LO in 
ChPT the ratios are given by 
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In practice, subleading terms in the chiral expansion 
are significant at the values of m and used in present 
simulations. Thus we have extended the calculation of 
the ratios to next-to-next-to-leading order (NNLO) in 
the power counting appropriate to the m ^ regime 
(usually called the "large cut-off effects" or LCE regime) . 
This power-counting is explained in Sec.|Tl] It differs from 
the usual continuum power-counting in that one-loop ef- 
fects are of NNLO, rather than next-to-leading order 
(NLO). The only NLO contributions are from analytic 
terms. 



The LO results are generalized in a fairly simple way. 
To describe this, we first define 'D{s,t,u) and S{s,t,u) 
as the infinite volume, PQ scattering amplitudes corre- 
sponding to the contractions of Figs.[lja) and (b) respec- 
tively [see Eqs. (|9| and ([lO| below] . We then observe that 
the coefficients of t/4M^L^ in Eqs. ^ and ^ are simply 
the LO values of V{4M^,0, 0) and S{4M^, 0, 0), the am- 
plitudes at threshold. What we show in Sec. |IV| is that 
the coefficients of t/iAM^L^), when evaluated to NNLO 
in the chiral expansion, continue to equal the infinite- 
volume threshold amplitudes, also evaluated to that or- 
der In effect, picking out the coefficient of the term 
linear in t in the ratios is, for the threshold amplitude, 
like performing LSZ reduction. 



Th e full NNLO results are given in Eqs. ( |104[ ) and 
( 105 1. We display here only simplified forms which show 



where Wq g are proportional to Wq g [see Eq. (26)], and 



f X f^. Thus these two LECs can be determmed from 
the terms linear in t. This requires that t be large enough 



^ This is true up to finite-volume corrections proportional to 
exp(— M,ri), which are generically present also in unquenched 
applications, and are usually small in actual simulations. 
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the essential features: 
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where V and S are evaluated to NNLO, and the results 
at threshold are given in Eqs. (106) and (107 1. In our 
"big O" notation, M"^/ p may also stand for any of the 
terms proportional to a^Wj^, as these are of the same 
order in the LCE regime. By "exp. suppr." we mean 
contributions which fall off exponentially with t, e.g. due 
to excited states. 

As can be seen from Eqs. ([S]) and ^ there are three 
expansions being used. First, there is the usual chiral ex- 
pansion (supplemented by powers of a) in the expressions 
for V and S. Second, there is an expansion in powers of 
t, as indicated in the last lines of ^ and ^ . Finally, for 
each power of t there is a sequence of subleading terms, 
as is shown in the middle line of each equation. 

Our calculation in Sec. |IV| yields the first non-trivial 
correction in each of these expansions, namely the chiral 
corrections to V and S and the contributions to the ratios 
proportional to / and t/L'^. The latter two are given 
explicitly in (104) and ( |105| ). We stress again that, if 



the numerators in the ratios were physical correlation 
functions, then, using the results of Ref. 'TSl, one would 
find that the t/L^ and t'^/L^ terms would be proportional 
to the square of the threshold scattering amplitude. This 
is not true in the PQ theory, but one can nevertheless 
calculate these terms. 

With the NNLO results (|3| and (|4| in hand, we can 
discuss our proposal for determining LECs in more detail. 
In order for the term linear in t to dominate over the 
quadratic term, it must satisfy 



t < pL^ 



{M^L){fLy 
M, 



(5) 



Since Mt^L ^ 1 and fL > 1 to avoid large finite- 
size effects, one sees that the constraint on t is rather 
weak. One may also try and fit the ratios including the 
quadratic term, and in this regard we note that the co- 
efficient of is given by a linear combination of M^/p 
and the LECs Wq and Wg, so that no new parameters 
are needed. 



Assuming that one can determine the coefficient of t, 
one must disentangle its dependence on and on 1/L 
in order to extract the LECs Wq and W^. Here it is help- 
ful that the 1/L correction itself depends on these same 
two LECs. At the least, this will allow an a posteriori 
estimate of the size of the 1/L correction. 

Finally, one must do a chiral extrapolation of the re- 
sulting threshold amplitudes, attempting to pick out the 
LO terms which give the desired LECs. From Eqs. (106) 
and (107), the general chiral behavior is, schematically, 

,3 



V,S 



a" + M2 4 
WM^{1 



a^{l 



logM,) + M^(l 



fogM^) 
-logM^). (6) 



The coefficients of the chiral logarithms are either fixed 
(for the continuum logarithm) or given in terms of all 
three W'^ (for the logarithms multiplying factors of a). 
The analytic terms involve many other LECs, however, 
both from the continuum theory and induced by dis- 
cretization errors. Thus it seems very unlikely that one 
will be able do more than a fit to the generic form given 
above and extract the LO and terms. This would 
allow a determination of W'^ and M^g, but not Wy. 

In light of this, we have devised an alternative method 
for determining W^, in which it contributes at tree- level. 
This requires studying a particular three-pion correlation 
function. Since this will be challenging to implement in 
a simulation, we describe the method only briefly in an 
appendix. 

We close this section by noting that other methods 
for determining Wg, Wlj and have recently been 
proposed. These use the eigenvalue distributions of 
the Hermitian Wilson-Dirac operator with Wilson-like 
fcrmions (which are sensitive to all three LECs) [101 E] ) 
the masses of pions and the scalar correlator in a 
mixed-action simulation with overlap valence quarks and 
twisted-mass sea quarks (which can determine W'^ and 
W^s) [IS] I s-^d the mass of the quark-connected neu- 
tral pion with twisted-mass quarks (which determines 
W^) [17]0 We think that pinning down the LECs will 
not be easy, and hope the method proposed here can 
contribute along with these other approaches. 

The remainder of this article is organized as follows. 
In the following section, we give a brief recapitulation of 
the pertinent details of PQWChPT, including the power- 



counting of the LCE regime. In Sec. Ill we present our 
results for the infinite-volume PQ scattering amplitudes, 
which we do for general momentum, and at NNLO in the 
LCE power counting. The description and calculation of 
the finite-volume correlations functions are presented in 
Sec. |IV[ which forms the core of the technical part of this 
paper. We include two appendices. The first contains de- 
tails concerning analytic NLO and NNLO contributions 
to the scattering amplitudes. The second describes our 
proposal for determining W^. 



It has also been shown in Refs. |10II11|[T7] that W'^ is necessarily 
negative. 



4 



II. PARTIALLY QUENCHED WILSON CHIRAL 
PERTURBATION THEORY 



In this section we define the required PQ scattering 
amphtudes and recah the essentials of WChPT for the 
partiaUy quenched theory. We consider a theory with 
two sea quarks, and introduce four valence quarks, and 
their corresponding ghosts, in order to define the desired 
amplitudes. All quarks and ghosts are degenerate with 
mass m. With this quark content, the chiral field 



E = exp {V2iTT/f 



(7) 



lies in the graded group 5f7(6|4). We work in the LCE 
regime where m ~ a^, so the Lagrangian is broken into 
leading and subleading parts as follows [12j : 



This is found by comparing quark-level Wick contrac- 
tions. We use this result to provide a partial check of our 
results for V and S by doing an independent computation 
of in SU{2) WChPT. 

We can also relate the PQ amplitudes to the general 
physical scattering amplitude. We recall that, for tt; -|- 
TTfe — )■ TT; + TTm, with i, k,l,m € {1, 2, 3}, one can write 



SikSlmA{s, t, u) + SaSkmA{t, S, u) + Sim5lkA{u, t, s) 



(13) 



The amplitude A is related to our amplitudes by 



LO 
NLO 
NNLO 



2 2 

p ,m,a 
p^a, ma, 

4 2 2 2 2 2 4 

p ,p m,ra ,p a ,ma ,a . 



(8) 



A{s, t, u) = 

V{t, s, u) ~ 5(s, t, u) + 5(w, t, s) + S{t, s, u) . (14) 



The partially quenched scattering amplitudes of inter- 
est correspond to the two contractions shown in Fig. [T] 
We label these, respectively, as V and S for double and 
single, referring to the number of loops appearing in 
quark-flow diagrams. The four valence flavors that we 
have introduced allow us to separate the contractions as 
illustrated in Fig. [2] The precise definitions are 

V{s, t, u){2tt)'^5^{p + k^p' -k')^ 

(7ri2(p')7r2l(-p)7r34(fc')'r43(-fc)) conn, amp (9) 

and 

5(s, t, u){2t:)^S\p + k-p' -k') = 

(7ri2(p')'^23(-p)7r34(fc')'^4i(-fc)) conn, amp i 

(10) 

where conn indicates that only pion-connected contribu- 
tions are included and amp indicates standard amputa- 
tion of external propagators. The subscripts on the pion 
field indicate their valence flavor and s, t and u are stan- 
dard Mandelstam variables, 



^-{p + k)\ t^-{p-p'f 



(11) 



with p, k, p' and k' all Euclidean four- vectors. 

Observe that the term "amplitude" only applies loosely 
here because of the unphysical partial quenching. T) and 
S do not satisfy unitarity constraints and only certain 
linear combinations give physical amplitudes. In partic- 
ular, the relation of PQ amplitudes to the amplitude for 
7r+ scattering is 



The result for A at NNLO in WChPT is obtained in 
Ref. |12j . and this relation allows us to compare our re- 
sults to those in that work. 

In order to calculate T) and S to NNLO we must in- 
clude all tree level and one loop diagrams generated by 
the LO Lagrangian (Clo) but only the tree level dia- 
grams from Cnlo and Cnnlo- The LO Lagrangian 

is laiiH] 



St) 



a^WjiJ:-!:'^)^ 



(15) 



Here angle brackets indicate a supertrace (strace), and 
X — 2Bom and a — 2Woa, with Bq and Wq leading- 
order LECs. The term linear in a has been removed by 
the standard redefinition of m [5 . We use the "small /" 
convention in which f-^ 93 MeV. 

When the field S is restricted to SU{2), the number 
of independent terms in £lo is reduced. In particular, 
the IFy term vanishes and the Wq and terms become 
proportional. As a result, the unquenched amplitudes de- 
fined in ( 12 1-( 14 ) can only depend on Wq and W^ through 
WL 



the combination 2Wq 
this, used in Ref. ITU 



A convenient definition for 



IS 



^7r+ {s, t, u) ~ 

V{s, t, u) + V{s, u, t) + S{s, i, u) + 5(s, it, t) . (12) 



W^ 

C2 = -32(2M^^ + I^^)-^ 



(16) 



The NLO and NNLO Lagrangians are [H [1^] 

-axW^6(S + St)2 

-axW^8(S2 + (St)2) 
+ a^— terms 



(17) 





(a) 



(b) 



(c) 



XX 



n 



(d) 



(e) 



(f) 



NNLO 



where 



■ xi5(9^sa^st(E + st)) 

X%(S + Et)2 

x2Lr(S-I]t)2 
x2L8(E2 + (st)2) 

a^7n— terms 
a^p^— terms 
a*— terms , 



OpQ = (9^E9,Eta^E5,Et) 



(18) 



(19) 



The a^, a^m, a^p^ and terms are discussed in ap- 
pendix |Aj 

As above, restriction to E e SU{2) reduces the number 
of terms in the Lagrangian. As a result, physical ampli- 
tudes can only depend on the following combinations of 
NLO and NNLO LECs, 



2Wi 



W^, 2We + Ws, 2Le+Ls, (20) 

2W^4 + M^5, 2^4 + ^5, (21) 

2Li-fL3, L2, (22) 

as well as on the LECs from a^, a^m, a^p^ and a"* terms. 
In particular, physical amplitudes must be independent 
of W^, Wy, L7 and Lpg because the associated La- 
grangian terms vanish when there are only two quarks. 



III. INFINITE VOLUME PARTIALLY 
QUENCHED SCATTERING AMPLITUDES 

In this section we calculate the on-shell PQ amplitudes 
through NNLO. Although, as noted earlier, the PQ the- 
ory is defined in Euclidean space, we can nevertheless 



(g) 



(h) 



FIG. 3. Classes of diagrams contributing to the PQ scatter- 
ing amplitudes through NNLO. Filled circles represent ver- 
tices from Clo while filled squares represent vertices from 
^NLO and /Innlo- 



analytically continue the amplitudes to Minkowski mo- 
menta and set them on shell. It turns out that it is these 
on-shell amplitudes which appear in the coefficient of t 
in the ratios of Fig. [2j 

The diagrams contributing through NNLO are shown 
in Fig. [sj At LO, only Fig. |3ja) contributes, and leads 
to the following results for the on-shell amplitudes: 



= 2w' 



1 
2/ 



2 (2*^0 



Here, Mq is the LO pion mass, given by 



X + f{2w', + w',), 



(23) 
(24) 

(25) 



where we have introduced rescaled, dimensionless LECs 

(26) 



Note that, because all quarks and ghosts are degenerate, 
all pions have the same mass, whatever their composi- 
tion. 



using the on-shell resu 



Combining Eqs. (23) and (24 1 according to (14 1 and 
t s + t + 



4Mo, we find 



MS 



P 



2w' 



Wo 



(27) 



This agrees with the result of Ref . [T2] . 

We now turn to higher-order contributions, consid- 
ering first the loop graphs of NNLO and then turning 
to the analytic contributions of NLO and NNLO. Non- 
analytic contributions to mass and wavefunction renor- 
malization arise from the "tadpole" diagrams exemplified 
by Fig. |3][b). Including the results from these diagrams 
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(and also the analytic analogs in Fig. [sj^h)), the form of 
the valence-valence propagator near the physical pole is 



tudes is 



(7rij(p)7rj,(-p)) 



(1 



1 



M2) 



(28) 



with ~ indicating that the two sides differ by terms reg- 
ular at the pole. This expression defines both the wave 
function renormalization bz and the physical mass Af^r. 
It is further convenient to define — Mq -|- 5M^. 
Evaluating the tadpole diagrams we find 



^^loop — 



2 
2 



Wo 



where 



Mi 



_Ml_ 
167r2/- 



is the standard tadpole integral. Here 



log ( A/2 /m' 



(29) 
(30) 

(31) 



V 



(2) 
loop 



(s,t,«) = -2<5z,„„pI?(") + ^/i 

+ (1/36) Ms) + hit) + hiu)] 
+ [-(10/3)u;^-2u;^ + 2u;'7]7i 
- (4/3XJ4(t) 

+ {l/3)w's[h{s)+h{u)-2h{t)] 
+ w'e^ [4l2is) + 4l2{u) + Uhit)] 
+ lOw'ew'shit) 

+ w'i Ms) +l2{u) + [17 /2)l2{t)] 



and 



(35) 



74— e. 



(2^: 



4-e 



(32) 



and the arrow indicates evaluation in a modified minimal 
subtraction scheme (MS) in which 'je — log(47r) — 1 is 
subtracted along with the pole 

The result for bMf^^^ is the same for all pions (as 
required by the graded symmetry group) including the 
physical, unquenched pions. It must thus contain dis- 
cretization errors proportional to the combination 2u'g + 
Wg, as we see is the case. The result (29) agrees with 
that given in Ref. |12| . 

Non-analytic contributions to the scattering ampli- 
tudes arise from the diagrams of Fig. [3f^c-f), from wave- 
function renormalization, and also from mass renormal- 
ization. The latter contribution is, at the order we work, 
only present in S , and arises kinematically from the LO 
result. To see this explicitly note that the general, off- 
shell form of the LO amplitude is 



c(2) 

loop 



{s,t,u) 



where 



3/2 
10s 



<5ACp-25z,<,„p5(")(s) 



13Af, 



18/2 

+ (1/18) [/7(i,u)+/7(u,t)] 

- [(ll/3)w^ + (3/2)^ + 2w'j] h 

+ (l/3)w;^ [2/4(5) - h{t) - h{u)] 

-il/3)w'sMt)+h{u)] 

+ Aw'ew'sMs)+l2{t) + l2{u)] 

-{5/2)w'iMt)+l2{u)] , 



(36) 



S^'J (s, t,u) = w', + -^ (2A/2 -2s + t + u) 



1 



6/ 



(33) 



When going on shell one sets p — k = p = k'~ to 
— A/2^ with A/tt the pion mass at the order being cal- 
culated. This imphes that s + t + u = 4(A/g2 + dlVP). 
Substituting this in (331 one finds a NNLO contribution 
to S of 



Ml 
3/2 



(37) 



S 



(2) 
(5M2 



25Af2 



(34) 



Combining this with the other loop contributions we find 
that the entire one-loop contribution to the PQ ampli- 



In each case, the results proportional to /i are from the 
tadpole diagram, Fig.[3]jc). 



The new integrals that appear, and their values after 
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subtraction and going on-shell, are 



his) 



^ [Fis) + 1 - log{M^/ti')] 
1 f -5s + t + u-2q^ -2q^ 



(38) 



(39) 



P 



2/i 



1 



h{s) = j^ J iSs+p^ + k'+q^ + q") 
jSs+p'^ +k'^ + q^ + q') 

(3s-4M|)2 10s - 16M2 



(40) 



P P 

{i{k' + qf - p'^ - k'^ - q^ - q^) 

3s2 + 3us/2 - l2M^s - 6M^u + 16M^ 
his) J, 



h 



4s + 3ii - lOMg i6M^ - s) i4M^ ~ s ~ 2u) 



P 



327r2/4 



where q = q + p + k. Thus aU integrals can be written in 
terms of Ii and a function defined in Ref. J2lI- 



Fis) = -as log , as 

CT„ — i 



(42) 



The analytic contributions of NLO and NNLO arise 
from Figs, [sjg) and (h). For the mass shift and wave- 
function renormalization we find 



SM^„ = -(l/2)(2w4 + W5)M^ + i2we + ws)x 
-il/2)i2U + ^5)M^ + i2^e + ^s)x 



(43) 



+ SAP 



ad 



6Zan = (1/2) [(2W4 + W^) + (2^4 + Cs)] + ^ad , (44) 

where 



Wk 



P 



6 



P 



(fc = 4-8) (45) 



are defined in analogy with Eq. ( 26 1 . SAp^^ and 5zad are 



the contributions from the "additional" terms, meaning 
a? , a^m, and a^p^ terms. They are given in Appendix 

El 

Incorporating the corrections from SAl^^ and 6zan we 
determine the NLO and NNLO analytic terms in the PQ 



amplitudes to have the explicit form 

Vgi^\s,t,u) = i2wi + W5)i-2w'e)+WiGit) 
+ '2x^6/ P 

+ 8LiG(t)2+4L2 (G(s)2 + G(u)2) 
- 4LpQ (G(s)2 + Git)^ + Giu)^) 
+ i2U + C5)i-2w'e)+UGit) 
+ 2x^6/P+Vadit), 
Si];'\s,t,u) = (2t.4 + W5)(-^ - Ap/p) 
+ W4s/i2p) 
+ 2xiwe + ws)/P 
+ 2L3 (G(t)2 + G(w)2) 
+ 4LpQ (G(s)2 + G(0' + G(u)2) 
+ i2^, + ^,)i-w's-Ap/P) 
+ UPi2p) 

+ '2xi^6+^8)/P+Sadis), 



(46) 



(47) 



where 



G(s) 



2M^ 



P 



(48) 



and Vadit) and Sadis) are given in Appendix [A| Note 
that, due to the great number of LECs, the analytic parts 
of the PQ amplitudes are poorly constrained. 



The final results are obtained by combining Eqs. (23 1, 



(35) and (46) for 



Vis, t, u) - + Vlllis, t, u) + ViY^ is, t, u), (49) 



and similarly combining Eqs. (24), (36) and (47) for 



Sis, t, u) = 5(0) is) + Sl'j(s, t, u) + Si'P^ is, t, u). (50) 



These may then be used in Eq. ( [T4| to find A. We find 
complete agreement with the result of Ref. ^2j . An inter- 
esting aspect of this comparison is that the terms linear 
in a in I? and S (which have the generic form aMg and 
as) cancel in A. This is only true when the leading o rde r 
amplitude is expressed in terms of Afp as in Eq (27). 



IV. FINITE VOLUME CORRELATION 
FUNCTIONS 

In this section we calculate, in PQWChPT, the corre- 
lation functions appearing in the ratios shown schemat- 
ically in Fig. [2] Specifically, we consider a cubic spatial 
box of length L with periodic boundary conditions, and 
assume that the time direction satisfies Lt ^ L and is 
effectively infinite. Since the PQ theory does not have a 
positive transfer matrix, correlation functions cannot be 
analyzed by inserting complete sets of states with positive 
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normj^ Furthermore, the scattering amplitudes of the 
infinite volume theory are not unitary. Because of these 
unphysical features, the standard relation due to Liischer 
between finite volume energy shifts and phase shifts [T3] 
does not hold — neither energy shifts nor phase shifts can 
be defined. Instead, if one wants to "measure" scattering 
amplitudes, one can compare results for Euclidean cor- 
relation functions (obtained using lattice methods) with 
the predictions of PQWChPT. Then the unphysical fea- 
tures of the underlying theory are matched by those of 
the effective theory. This strategy was introduced in 
Ref. |15j to study pion scattering in the quenched ap- 
proximation. 

To approach infinite-volume scattering as closely as 
possible, we consider external pion fields with definite 
three-momentum, leaving only their time coordinates un- 
transformed. We additionally restrict ourselves to the 
simplest case of pions at rest. Specifically, we calculate 
the following four-point correlators: 

Cvit) = (^12W^34(<)^=2l(0)i=43(0)> (51) 
Csit) = (^I2(0^34(t)i=23(0)i=4l(0)) , (52) 

with 



7?(i) = / d^xT:{x,t) . (53) 

The subscript on the integral indicates integration over 
finite volume. These correlation functions are, roughly 
speaking, the finite- volume analogs of the infinite- volume 
unamputated scattering amplitudes. In order to more di- 
rectly access these amplitudes, we take the ratio of these 
correlators to the square of a single-pion correlator, which 
is, roughly speaking, the analog of amputation; 



C^{t) = (i=12(t)^^2l(0)) 



(54) 
(55) 



These are the ratios shown schematically in Fig. [2j In 
this note we will always take t to be positive. 

It is instructive to make contact with the more familiar 
resu lts i n a physical (i.e. unquenched) theory. Adapting 
Eq. ( 12 1 to finite volume, we find that the 7r+7r+ correla- 



tion function is related to our PQ correlators in a simple 
way: 



C.+ (i) = (^+(i)^+(i)^" (0)^^(0)) 
= 2Cv{t) + 2Cs{t). 



(56) 



It follows that 
Rv{t) + Rs{t) = Ze"(^-+-+"2*^'''* + exp. supp., (57) 



where 



is the energy of the lightest 7r+7r+ state 



with zero total spatial momentum, and the exponentially 



As seen explicitly in the transfer matrix obtained in Ref. |21| 



suppressed terms come from states with higher energies. 
At LO in WChPT one finds (as we will show below) that 
the overlap factor Z is unity, so that 



Rv{t) + Rs{t) Kl~SEt 



6E 



2AU 



(58) 



This approximation is valid as long as t is large enough 
that the exponentially suppressed terms are negligible 
but also small enough that the linear term dominates 
the Taylor expansion of the leading exponential. 

Liischer's result relates this shift to the infinite volume 
TT"*" scattering amplitude at threshold [13 : 



SE = -A 



th 



1 



0(e-*^"^) , 



th 



1 



IGttA'UL 



0(1/2.2 



(59) 



where ci = —2.837297 is a numerical constant. The form 
of the l/L^ correction, and some higher order terms, are 
known DRT , but we do not show them as we will not con- 
trol the corresponding terms in our calculation. Note also 
the presence of exponentially suppressed finite-volume 
corrections. 

Although the PQ ratios Rv{t) and Rs{t) do not behave 
as a sum of exponentials, what we can take over from 
the analysis of R„+{t) is that it is useful to determine 
the coefficient of the term linear in t. In the remainder 
of this section we determine the PQWChPT prediction 
for the linear terms in Rv{t) and Rs{t)- We work to 
NNLO in the momentum power counting of Eq. Q and 
do so controlling not only the leading terms but 

also the corrections. We can also control a subset 

of the finite volume corrections proportional to e"*^"^, 
but will not do so systematically. 

The PQWChPT diagrams which contribute to the or- 
der we work are shown in Fig. [4] These are the same 
diagrams as for the infinite-volume amplitudes. Fig. [3j 
except for the addition of the disconnected diagrams (a- 
c). Our description of the calculation is broken into 
three subsections: leading-order results, analytic NLO 
and NNLO contributions, and NNLO results from loop 
diagrams. 



A. Leading-order Results 

The pion-disconnected diagram of Fig.Qa) contributes 
only to Cx){t), with the result 



Op 



'{t) = C^{t) 



Mot 



2Mo 



and thus 



R. 



(disc) 

V 



it) 



1, Rf''\t) 



0. 



(60) 
(61) 

(62) 



The factor of in (^^(i) arises because both ends of the 
propagators are integrated over space. Note that, within 
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Q- 



(a) 



(b) 



(c) 




(d) 



(e) 



(f) 



(g) 



XX 



(h) 



(i) 



0) 



FIG. 4. Classes of diagrams contributing to the correlators 
C-oit) and Cs{t) in ChPT. Diagrams in which initial or final 
pions are interchanged are not shown separately. Notation for 
vertices as in Fig. [3] 



WChPT, there are no explicit contributions from excited 
pions at any order. The effect of these states appears 
through contact terms proportional to S{t), which arise 
first at NNLO. Since we always consider t large enough 
to remove terms which are exponentially suppressed, we 
need not include such contact terms. 

The tree-level, pion-connected diagram, Fig.|4]jd), con- 
tributes to both Ct) and Cs- Because the four pion 
sources have p = 0, only pions at rest contribute. Thus 
the pion propagators attaching to the vertex are either 



g-Mo|ti| g-M(,|t-ti| 

or 



2A/n 



2Mn 



(63) 



where ti is the time of the vertex. These propagators 
have no factors of since only one end is integrated over 
space. There is, however, a factor of from integrating 
the position of the vertex over space. 

The correlator Cv{t) is the simplest to consider, be- 
cause the LO vertex comes only from the Wq term in 
£lo [see Eq. (15l] and therefore contains no derivatives. 
Based on the discussion above, we find (recalling that t 
is positive) 



a 



(conn,0) 
V 



-2Mo|ti| -2Ai?o|t-ti| 

it) = 2w',L' / dt, — — (64) 



g-2Mot 



(2Mo) 



+ t 



2Mo 



(65) 



The factor of t arises from the region < ii < t in 
which the vertex lies between the two sources. The t- 
independent terms come from ti < and ti > t, where 
the contribution drops exponentially. 

At this point, it is useful to incorporate a result from 
the calculation of subleading orders. The diagrams which 
correct pion propagators, Figs.[4][b), (c), (e) and (f), form 



part of the geometric series which changes 



,-Mo\t\ 
2Mn 



to 



,-M^\t\ 

~2M7 



(66) 



Here is the physical pion mass to the order we are 
working, and = 1 — 5z is the wavefunction renormal- 
ization. For the moment we incorporate only the mass- 
shift, returning to the effect of ^ 1 below. Using the 
new propagators in both numerator and denominator of 
the ratio R-pit), we find 



R 



{conn,0) 



it) 



1 



/ 1 



1 f 1 



+ t]2w'. 



6 J 



4M|L3 \2AU 



(67) 
(68) 



The utility of the ratio is that the overall exponentials 
cancel — as noted above, this corresponds to amputation 
in a continuum calculation. The physical interpretation 
of the t term is that the pions can interact at any inter- 
mediate time, while the suppression arises because 
the zero-momentum pions must overlap in order to in- 
teract. As shown in the second line, the coefficient of 
t is, aside from the kinematical factor 1/(4M^L'^), the 
LO PQ scattering amplitude. This, together with the 
disconnected contribution from ( [62| ), gives the result ^ 
quoted in the introduction. 

A similar analysis holds for Rs{t), except that we must 
now deal with the momentum dependence arising from 
the kinetic term in the LO Lagrangian (15). If one eval- 



uates the diagram in position space, the derivatives in 
the vertex act on the pion propagators of Eq. ([63]) (with 
Mq — >■ M^). Thus only time derivatives contribute, and 
they yield zLM.^. Consider first < ti < t, i.e. the ver- 
tex lying between the sources. Derivatives acting on pion 
propagators originating at times and t then give —M^^ 
and -t-Afjr, respectively. This implies that s = 4M^ and 
t = u = 0, i.e. on-shell kinematics at threshold. For 
ti < {ti > 0), by contrast, all derivatives give +AlTr 
(— M^), and one obtains, in both cases, the amplitude at 
off-shell kinematics: s = t = u = 4M^. The final result 
is 



R 



(conn,0) 



1 



(m-A/oV3/^) 



4M2L3 



Wo 



2M^ 



3/2 



(69) 



1 



4Af2L3 



2Af, 



5(°)(4M2,4Af2,4Af^) 



-5(")(4Af2)t 



(70) 



This result, together with the vanishing of the discon- 
nected contribution, is reported in Eq. ([2]) of the intro- 
duction. Note that, since in ^ we are quoting a LO 
result, we can set Mq = Af^. 



10 



As a check on our results we consider the sum i?p (t) + 
Rs{t). According to the discussion above, this should be 
« (1 - SEt), with 6E ^ ~A^^+/{4:M^L^) [the LO term 
in Eq. ([59])]. Wc find (setting'^Afo = M^) 



R 



l + O 



\t) + R: 
1 

Z3 



(conn,0) / , \ 



1 



2wr 



Mi 
P 



(71) 



The coefficient of t/ [AM'^L^) is indeed the 7r+ scattering 
amplitude at threshold. The corrections to the t- 

indcpendent terms (which are not shown in detail) are 
the first correction to the Z-factor for the two pion state. 



B. Analytic NLO and NNLO contributions 

Analytic NLO and NNLO contributions arise from 
Figs. |4][b), (e) and (g). Diagrams (b) and (e) give the 
analytic parts of mass and wavefunction renormalization, 
contributions which are identical to those in infinite vol- 
ume. The effect of mass renormalization has already been 
discussed above. Wavefunction renormalization partly 
cancels in the ratios, leaving a factor of = (1 — 2Jz), 
exactly what is needed to renormalize the amplitude as 
in infinite volume. 

The analytic contributions to the vertex, Fig. [4]Jg) , can 
be analyzed by a straightforward generalization of the 
method used for the momentum-dependent contribution 
to the LO vertex in Rg- We find 

t 



R. 



(1.2) 
"D ,an 



it) 



const 



R 



const 



t 



25(1/) (4M2, 0,0), (72) 
5(J,'2)(4M2,0,0), (73) 



where "D^an^^ and iSan^' are the infinite-volume analytic 



contributions to the two amplitudes given in Eqs. (46) 
and (47 1. We note that, at threshold, these amplitudes 
are just linear combinations of a"^, aMj, 



and 

with independent coefficients. The results (|72|) and 



( 73 1 are the natural generalizations of Eqs. ( 68 ) and ( 70 ) 



Note that we only keep track of the terms linear in t, since 
these are proportional to the desired PQ amplitudes at 
threshold. The constant terms involve off-shell ampli- 
tudes. 



C. NNLO results from loop diagrams 

In this section we extend the calculation to the one 
loop diagrams which appear at NNLO, focusing on the 
coefficient of the terms linear and quadratic in t in the 
ratios RT>,s{t)- At one loop there are 3 types of contribu- 
tions, which we discuss in order of increasing complexity. 
First, there are tadpole diagrams, shown in Fig. |4]jc), (f) 
and (h). Second there are t and u-channel loops, exem- 
plified by Fig.|4|^i). And, third, there are s-channel loops, 
as shown in Fig. |4|^j). 




d 



(g) 



(h) 



(i) 



FIG. 5. Classes of diagrams contributing to the correlators 
C-vit) and Cs{t), arising from tt^ and tt^ terms in the external 
operators. All vertices are from the LO Lagrangian. 



1. Tadpole diagrams 

Tadpole diagrams on external legs [Figs. |4|^c) and (f )] 
renormalize the external pion propagators, contributing 
to dlVP and Sz as described for the analytic terms. In 
this case, however, there is a difference compared to in- 
finite volume, namely that in the tadpole integral Ii one 
should use the finite volume pion propagator. This gives 
rise to corrections which are suppressed by powers of 
exp(— A/jr^), as can be seen by implementing the peri- 
odic boundary conditions using images. We assume such 
corrections are negligible. 

Tadpole diagrams attached to the vertex [Fig. Qh)] 
are also simple to incorporate. As long as < ii < i, 
they multiply the tree-level vertex by a factor that is 
independent of the vertex position, leading to the usual 
term linear in t. The coefficient of t is proportional to the 
tadpole contributions to the threshold amplitudes, and is 
exactly that needed to maintain the forms of Eqs. (72) 



and ( 73 ) , with D and S now including one-loop tadpole 
contributions. Again, there are exponentially suppressed 
volume corrections which we assume negligible. 

At this stage, it is appropriate to mention that there 
can also be tadpoles arising on the external operators at 
times and t. This is because, in practice, if one uses 
a local pseudoscalar operator at the quark level, then 
it maps into the chiral theory at LO as qklbQji^) ~^ 
c(S — S^)jfc(x). The constant c is known but unimportant 
here, since it cancels in the ratios. This chiral operator 
expands to a term proportional to ■Kjk — the operator we 
have been using — but with corrections proportional to 
{T^^)jk/P and {T^^)jk/f^- These corrections give rise, 
at the order we are working, to the diagrams of Fig. [5] 
As we now explain, however, all of these diagrams lead 
to contributions subleading compared to those we are 
keeping. 

Tadpoles associated with external operators [Figs, [s]^ a) 
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and (b)] cancel in the ratios. The loop diagram of 
Fig. [sjjc) does not give rise to a term linear in t because 
the four propagators to the left of the vertex lead to a fall 
off of at least e"^*^"*! . The tree-level diagram of Fig.jsjjd) 
also does not contribute to the coefficient of t — it gives 
a contribution to the constant term. The same is 

true for the corresponding tadpole diagram, Fig. [Sjje). 
The only diagram in this class that does give a contribu- 
tion linear in t is Fig. [sjjf). This arises when both pions 
in the loop have q = 0. This contribution is, however, 
of size t/L^ in the ratiosj^and thus is sublcading to the 
1 /i* terms that we are aiming to control. 

There are also analytic corrections to the external op- 



erators, e.g. terms proportional to d'^TTjk and 



These lead to corrections which cancel in the ratios. 



2. t and u-channel loops 

The t and w-channel loop diagrams have the form of 
Fig. Qi). Note that, since both external pions at time t 
are summed over all space, there is no difference between 
the t and it-channel loops. In order to avoid confusion 
between the two uses of t we will couch our discussion in 
terms of the u-channel. 

Figure gives rise to a contribution proportional to 
the time separation t as follows. Although there are two 
vertices (at times ti and 12), when they are pulled apart 
in Euclidean space there is an exponential suppression, so 
the dominant contributions occur when |ii — < ^/M^. 
The loop collapses to an effective vertex, which, when 
integrated over the intermediate time, leads to a factor 
of t. If either of the vertices is outside of the region 
< ti,2 < t, it is easy to see that one does not get a term 
linear in t. 

We now turn these words into a concrete evaluation. 
We consider first a contribution in which both vertices 
have no derivatives. In infinite volume, this leads to those 
terms in Eqs. (35 1 and (36) which contain the integral 
/2(m). For definiteness, we consider the contribution to 



■^[oop '^^ 4?«g^/2(M). The corresponding contribution to 
the finite-volume correlator is 



Cvit) D 



1 



(74) 



h{t)^Y^ dh dt2 / G^{xi,X2r. (75) 

^ Jo Jo Jxi Jx2 



^ The occurs because the numerator in the ratios is inde- 

pendent of L, while the denominator is proportional to L®. The 
numerator is independent of L because the three external propa- 
gators have one leg summed and are thus L-independent, leaving 
two point-to-point propagators in the loop (each proportional to 
and two vertices integrated over space (each giving a factor 
of L3). 



Here G^(xi,X2) is the Euclidean pion propagator, which 
is related to 0-^(1) by 

C^{ti-t2)= [ [ G„{xi,ti;x2,t2). (76) 

J Xi J X2 

The contribution to the ratio is 



Rv{t) D 



1 



4</2(t) 



(77) 



It is straightforward to evaluate l2{t) explicitly, and one 
finds that it consists of a term linear in t up to corrections 
falling as e"^*^"*. A simple way to pick out the coefficient 
of t is to take a time derivative and then send t — ^ 00. In 
this way we arrive at 



Rv{t) D 



t 



^w'il'2icx^), 



(78) 



for sufficiently large t. 

To evaluate /2(^) note that l2(t) can be rewritten 

as 

l2{t)^2f dhf'dt2f G^{QM;x2Mf, (79) 

Jo Jo Jx2 

where we have used translation invariance and the sym- 
metry G^(x,0) = G7r(0, x). Thus its derivative is 



/^(t) = 2 / dt2 / G,(0,0;X2,t; 



(80) 



Given the exponential fall-off of G^r, this integral asymp- 
totes to its large t value once t ^ l/M-^. The factor of 2 
can be traded for an extension of the integral to negative 
values of t. In this way we find 

4(c»)= f dt2 ( G,(0,0;f2,i2)' (81) 

' E '5.(g)^ (82) 



27r L3 



q—27Tn/ L 



where Gjr is pion propagator in momentum space. This is 
simply the finite volume version of l2{u = 0) = J^l/{q'^ + 



M2)2: 



/^(oo) -/2(u = 0) + O(e-*^'^) 



(83) 



where images can be used to see that the finite-volume 
corrections fall exponentially. 

The conclusion of this analysis is that, for sufficiently 
large t and L, the diagram which leads to a contribution 
to I'joop of 4?x;g^/2(u) contributes to the finite- volume ra- 
tio as 



Rvit) D 



t 



(84) 



Thus, once again, the coefficient of t/ (AM^L^) is the 
infinite-volume scattering amplitude at threshold. The 
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same argument goes through identically for I2 (u) contri- 
butions to iSj^jp, and for contributions to both amplitudes 
proportional to l2{t) {t here the Mandelstam variable). 

A similar argument also holds for the contributions to 
'^loop ^'^'^ ^loop proportional to the integrals l4{t), ^4(u), 
leit), leiu), lT(t,s,u) and I'j{u,s,t). These arise when 
one or both of the vertices are from the kinetic term in the 
Lagrangian. By similar manipulations to those above, 
one finds that the contribution proportional to t involves 
exactly the integrands of the infinite volume forms ( 39 
41), evaluated at threshold, but the spatial momentum 
integrals are again replaced by finite-volume sums. The 
short-distance divergences are unaffected by the finite- 
ness of the volume, so the regularization and subtrac- 
tions are unchanged. Thus the manipulations relating 
these integrals to I2 and Ii go through, as always up to 
exponentially small volume corrections. The net effect is 
that the forms of Eqs. ( 72 ) and ( 73 ) are maintained, with 



T) and S now including t and u-channel contributions. 

There is one subtlety for terms with two or more 
derivatives acting on the same internal propagator. To 
illustrate this point we define the integral 



1 



t pt 
dti I dt2 



J X\ J Xi 

2 



'x{diGT,(xx;x2)\G^(xi;x2)- (85) 



Although the integrand of ( 85 1 is an even function of x\ 



X2 , one must be careful in carrying out the manipulations 
which give the analog of Eq. ( 79 ) . The issue is that 



has a cusp at t\ — t2 and therefore the time derivatives 
give a delta-function: &(t\ — t-i). Carefully including the 
region about the delta-function, we find 



Iqi it) — lim / di\ 



2 / dt2 





tl+e 



dU 



ti-e 



X [dl^G^{d,h;x2,t2)]G^{6,h;x2,t2). (86) 

We now proceed as above, taking the time derivative and 
sending t 00 to deduce 



/a2(oo) = lim 



dt2 + / dt 



dU 



X / [dlfi^{Q-x2)]G^{Q;x2) 



(87) 



d''x2[dlG^{0;x2)]G^{Q;x2). 



As above, we are left with a single integral over all of 
space-time (with space finite). Note however that if the 
middle integral were absent then the final equality would 
not hold. Proceeding from Eq. (88), it is straightforward 



to see that 7^2 (00) is equal to the corresponding finite- 
volume amplitude at threshold. 







ti t2 t 

(a) 



t2 ti 

(b) 



t, 



ti 



t u 



(c) 



(d) 



FIG. 6. Time-orderings of the s-channel loop which give rise 
to contributions to R-D,s{t) linear in t at large times. The 
time-ordering (a) also give rise to a quadratic term, as dis- 
cussed in the text. 



3. s-channel loop 

The s-channel loop diagram is shown in Fig.Blj). This 
diagram leads to the terms proportional to l2{s), ^4(5) 
and Ie{s) in the infinite volume amplitudes. 

We begin by analyzing the case in which neither vertex 
has momentum dependence, which leads to the integral 
12(5) in infinite volume. For definitenes s we focus on the 

w'ihis) contribution to V^l^ in Eq. (35). The corre- 
sponding contribution to Cx>{t) is 



Cv{t) D 

~'2 



g-2M,|ti| g-2M,|t-t2| 

^2^[G.(.r,.2)]^^,^. (89) 



There are many choices of the ordering of the four times 
0, ii, t2 and t, but only the four shown in Fig. |6] lead 
to terms linear in t. The others give constants or terms 
which fall exponentially with t. 

In fact, the time ordering of Fig. [6][a) also gives a term 
quadratic in t. This occurs when both pious in the loop 
are at rest. The integrand is then independent of ti and 
^2, so that the contribution to the ratio is 



(90) 



In a physical theory, this is one of the contributions 
which builds up the quadratic term in the expansion of 
exp(— (5i?i) in Eq. (57), and we will use this below as a 
check on our result. First, however, we determine the 
quadratic terms arising from the finite volume versions 
of the integrals /4(s) and /6('S). 

These integrals have one or both vertices from the ki- 
netic term in Llo (Eq. ( 15 )). For our particular kinemat- 
ics the derivatives are simple to evaluate, since all pious 
are on-shell and at rest. The derivatives in the kinetic 
vertices thus give exactly the results that we obtain in 
infinite volume when working at threshold and on-shell 
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(s = AMI, t = ' 
substitutions, we can use Eqs. (39 



0). This means that, with these 



and ( 40 ) to relate the 



finite volume versions of I^is) and Iq{s) to that of l2is)- 
The Ii terms are to be dropped, since they arise from the 
q^O part of the loop. 

Putting this together, we find the total quadratic terms 
to be: 



Rvit) D 



AMI 
3/2 

_ AMI 
' 3/2 



- Aw, 



12 



1 



4M2L3 



\AM^L^ 



^ (91) 



(92) 



Thus, if one could determine the coefHcients of the 
quadratic terms in a simulation, one would gain addi- 
tional information concerning w'q and Wg. Note that the 
coefficients of the quadratic terms are not the squares of 
those of the linear terms give in Eqs. (67 1 and (69 1. This 



is another indication that the PQ theory is unphysical. 
The sum of the ratios is, however, physical and should 
be an exponential. Indeed, we find that the quadratic 
term in Rx) (t) + Rs {t) is exactly that needed to give the 
quadratic term in the expansion of Eq. (71). This pro- 
vides a non-trivial check on our results. 

We now return to the terms linear in t. We can first 
dispense with the time-orderings of Figs, ^c) and (d). 
Here a linear term arises only when the pions in the loop 
have q = 0. This means that the contribution to Cx},s{t) 
is of order and thus that to Rv,s{t) is order 1/L^, 
which is higher order than we are controlling. The same 
suppression holds for the contribution of pions with q = Q 
in time ordering of Fig. [6]Jb) . Thus we are left to consider 
the (f 7^ contributions to the time-orderings of Figs.[6]^a) 
and (b). 

Using Eq. (89), the time ordering of Fig.[6]ja) gives the 
following contribution to R'p{t): 



Rvit) D w'i 



L6 



dU 







dti 







-2Af„(ti-t2) 



(2M,)2 



(93) 



where the bar over the pion propagator indicates that 
the q = mode has been removed. We are interested in 
the coefficient of t/{AM^L^) for large t, and so we again 
take a derivative with respect to t and send i — > oo. The 
result, after some manipulations, is w'^Ia, with 



dti 



(94) 



Note that the growth of the exponential as ti becomes 

2 

more negative is overwhelmed by the decrease in 

A similar calculation for the time-ordering of Fig. [6[b) 
yields, for the coefficient of t/{AM^L'^), the result w^h 
with 



h = 



dh 



-2M, 



''G^ixi,Of 



(95) 



Here both the exponential and the G^ factors decrease 
as ti becomes larger. Combining the two time-orderings 



we find w'f times 



d^x, e-^^'-''G^ixi,Oy 



his = AM^) 



2\FV 



(96) 
(97) 



As noted in the second line, this integral is the finite vol- 
ume version of the corresponding infinite volume integral 
at threshold, expressed in position space. The factor of 
exp(— 2M7rti) simply leads to the injection of the physical 
threshold four-momentum through the loop. 
If we could ignore the difference 



Sh 



AM. 



2\FV 



Us = AMI)., 



(98) 



the result (971 tells us that, in the coefficient of 
t/{AMlL^), all the NNLO corrections to V and S pro- 
portional to l2{s) appear exactly as in infinite volume. 
As we will see shortly, however, 5I2 oc 1/L (plus expo- 
nentially suppressed terms) so we do need to calculate 
it. 

First, however, we extend the calculation of the coeffi- 
cient of t to s-channel loops with momentum-dependent 
vertices — the loops which give rise to the integrals /4(s) 
and I(i{s)- The steps outlined above lead to the same 
combined integral ( 96 1 , except that there are two or four 



derivatives acting on the various factors in the integrand. 
As was the case with the it-channel diagrams, one must 
carefully handle the discontinuity at ti = ^2 when two 
derivatives act on an internal propagator. The end re- 
sult, however, is still as claimed. This allows us to rewrite 
all contributions to the time dependent correlator which 
are linear in t as the corresponding contributions to the 
finite volume threshold amplitudes. Next we may relate 
these to the l2{s = AM^)^^ and /f^ using the same ex- 
pressions, (39) and (40), as hold in infinite volume (with 
s = AM^). From this follows that the difference between 
the finite and infinite volume forms is exponentially sup- 
pressed. Thus, aside from the SI2 terms, we find that the 
coefficient of t/ (AM'^L^) generated by s-channcl loop di- 
agrams is just full infinite-volume contribution from the 
same diagrams evaluated at threshold. 

Our final task is to evaluate 6l2- Standard manipula- 
tions lead to the following expression 



Sh = 



1 

Z3 



d\ 

(27r)3 



1 



(99) 



where -B^ 



VW] 



2 -I- M2 . Since the U V divergences can- 
cel in the difference between sum and integral, we can 
introduce a regulator term exp(— a[g]2), as long as we 
send a — )■ O"*". With this regulator in place, we can use 
Liischer's summation formula [I'd] (in the particular form 
quoted in Ref. ^T5\ ) 



L3 



d\ 

(27r)3 



fif) ci/(0) r(0) 
At:L L3 



q" 



(100) 
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where C\ is the constant given after Eq. ( 59 ) . This result 



is vahd up to exponentiaUy small corrections as long as / 
and all its partial derivatives are square integrable (as is 
the case with our regularized sum). Applying this result, 
we find 



ci 



1 



167rM^L 8(A/^L)3 



(101) 



At the order we are working, we need keep only the 1/L 
term. 

Collecting all the 1/L terms, we find their contribution 
to the ratios to be: 



AMI 
3/2 



4w 



1 



Clt 



(102) 

1 Clt 



AMi 



3/2 J AM^L^ m-KlVUL 



(103) 



The combinations of LECs here are the same as in the t"^ 



terms, Eqs. (91) and (92 1, because the both arise from 



the /2(s) integral in finite volume. 

A check on this result is that the t/L'^ term in the 
7r"'"7r+ correlator Rxi{t)+Rs{t) agrees with that obtained 
with Liischer's general formula (59). 



Here D is the full NNLO amplitude. The corresponding 
result for Rs{t) is 



Rs{t) = O 



Ml 1 



^ AMI 



5(4M2,0,0) 



t 



AMIL^ 



3/2 



1 



Clt 



1 



^2' 



4A/2L3 iQt:M^L yAM^L^ J 2 
''Ml 1 

X {l + O 



P M^L 



O 



Ml t 
P M^L^ 



+ 0(6"*^"-^) + exp. suppr. 



(105) 



Note that the last three lines of Eqs. ( 104) and ( 105 ) are 
identical. 

Thus, if one can measure the coefficients of i/i^, one 
obtains the corresponding infinite volume PQ threshold 
scattering amplitudes at NNLO. In fact, it is plausible 
that this holds to all orders, since the analysis above 
shows how picking out the t term corresponds to LSZ 
reduction in infinite volume. Of course, it is non-trivial 
to pick out these coefficients, but our results provide the 
coefficients of the leading competing terms — those pro- 
portional to t/L^ and t^/L^. 

For completeness we give the PQ results for the thresh- 
old scattering amplitudes. Using the results F(4M2) = 
0, F(0) = -2 and the notation Q = log{Ml / fi"^) / {IGtt'^) 
these are 



D. Summary of results 



Collecting the results from this and the previous sec- 
tion, we find 



Rv{t) = 1 + 



Ml 1 



P M^L^ J AM^L^ 



V{AMl 0,0) 



AMI 



3/2 



Aw', 



Clt 



( 1 



^i2' 



I 4M2L3 i67rM^L ^4^22.3 ) 2 

(Ml 1 \ 
V P M^LJ 



O 



Ml t 



n 3N 



P MIL^ 



+ 0{e ^''''^) + exp. suppr. , 



(104) 



V{AMl, 0, 0) = 2w'e + Vili^'> {AMl,0, 0) 



1 

167r2 



M: 



m: 



,f2 



2P / 



Q 



^ M 4 M2 

|^ + ^(4<-K2u;^-5z.^) 



22m;^ 



10w>g - (21/2)w^2 



(106) 



and 



/\/r2 r 71 ,r2 

5(4M2, 0, 0) = ^ - + Silf\AMl,0, 0)--^ 



1 



167r2 



A,f 4 9 /\,f 2 



Q 



- + 7^(2^6- 2^7 + 8u;^) 
- 12wgWg -1- 5wf 



(107) 



We do not give explicit expressions for I?i^^''(4M2, 0, 0) 

and 5i^,'^^(4M2,0,0) as they can be readily determined 
from Eqs. (liel) and (IItI). 
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Appendix A: Additional contributions to T> and S 

In this appendix we analyze the additional contribu- 
tions to the PQ amplitudes that come from the a'^ terms 
in the NLO Lagrangian and from the ma^, and p^a^ 
terms in the NNLO. We begin by enumerating all of the 

terms allowed by chiral symmetry: 

/:„3 - a^(E + I]n3 



(Al) 



f (I]t)3) 



We use ^ throughout this appendix to indicate that the 
two sides are equal with an independent LEG multiply- 
ing each term. The key observation is that these terms, 
when expanded in tt, produce (tt^), (tt^)^ and (tt"') with 
independent coefRcients. The particular forms of these 
coefficients in terms of the unknown LECs provides no 
useful information. 

Both the ma? and a* sectors generate the same three 
pion terms with independent coefficients. The argument 
for the ma? terms is identical to that for the a^ sector, 
since the spurious A and x transform in the same way. 
For the a* sector we can show this result by display- 
ing three chiral operators which give linearly independent 
contributions to the three pionic operators. An example 
is 



aHY. 



(A2) 



where the dots indicate additional terms which we do not 
need to enumerate. 

It remains to consider the p^a^ terms. Starting on 
the level of the pion fields, we first list all two-derivative 
quadratic and quartic terms. This is done without re- 
gard to chiral symmetry, using only that (tt) = 0. The 
resulting set is 

(5p7r7r)(9p7r7r) 

(a^7ra^^^2^ (A3) 

(9^7r TT d^TT tt) 
(a^7ra^7r)(7r2) . 



Because this is the maximal set of pionic terms to the 
order we are working, it is sufficient to show that thep^a^ 
terms allowed by chiral symmetry produce the entire set, 
with independent coefficients. This is indeed the 
follows from 



+ a2(9^I]a^S(Et)2) + ... . 



(A4) 



These terms are enough to independently give the set 
®- 

From these results, we can determine the contribution 
of all a^, a^TO, a'^ and p^a^ terms to the mass and wave 
function renormalizations and to the PQ amplitudes. We 
find 



SZad 

Vad{t) 
Sad{s) 



(?Ml ^'a?^ ^a" 

c?t + o^Mq +a^ + a^ 



c?s + c?Ml + + a^~ 2SzadS'^°^ (s) 

r2 //oj-2\ 



(A5) 
(A6) 
(A7) 

(A8) 



Appendix B: Determining Wj 

In this appendix we sketch a method for determining 
the LEG in which its contribution appears at tree- 
level. The method is by no means unique, but it is the 
simplest approach we have found within the context of 
pion scattering. 

Expanding out the term in the LO chiral La- 
grangian (15 1, the first non- vanishing term is 



Wj 

I8p 



-3\2 



(Bl) 



To get a tree-level contribution from this vertex one needs 
six external pious. We propose calculating the following 
finite-volume correlation function 

Cs^t) = (^=12 (0)^23 (0)^45 (0)^31 (i)^56(i)^64(0) (B2) 

and then forming the ratio 



a(t)3 



(B3) 



Here iTjk are the p — fields defined in Eq. ( (53|) a nd the 
single-pion correlator GT^{t) is defined in Eq. ( |55p . The 
subscripts on the pion fields indicate valence flavors, of 
which there must be six. We consider only i > in the 
following. 



The choice of valence fields in ( B2 ) allows only a sin- 
gle quark-level contraction, shown in Fig. [7j By con- 
struction, this contraction has two quark loops, in order 
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t=0 



t 





(a) 



(b) 





(c) 



(d) 



FIG. 7. Quark contraction for the correlator Cs-n of Eq. (|B2 



to match with the double-strace six-pion vertex of (Bl ). 



This wiU be a more challenging contraction to calculate 
in numerical simulations than those of Fig. [l] because 
of the "source to source" propagators at times and t. 
Nevertheless, with recent advances in calculations of "all- 
to-all" propagators, we expect that the calculation should 
be feasible. 

We have written the correlator in terms of the pion 
fields from the chiral Lagrangian, but, as discussed in 
Sec. |IV C H in practice one would use a quark- level pseu- 
doscalar field, qklb<lj- The corresponding chiral operator 



(E - Et),, = c 



1 

3/ 



(B4) 



where the constant c is known but not needed. It turns 
out that in this calculation, unlike that in the main text, 
the TT'^ part of the interpolating operator contributes at 
leading order to the quantities of interest. This means 
that it is essential for the present method to use (a dis- 
cretization of) a local pseudoscalar bilinear to create the 
pion fields, and not, for example, a non-local operator. 

Our choice of flavor indices significantly restricts the 
diagrams that can appear and the vertices that con- 
tribute. There is no pion-disconnected diagram, and the 
two tree-level diagrams which contribute are shown in 
Figs. |8]^a) and (b). Note that the four-pion vertices in 
(b) must be attached to the external legs as shown, other 
possibilities being forbidden by the flavor indices. There 
are also contributions involving the tt'^ and tt^ parts of in- 
terpolating fields, Eq. (B4|. For reasons explained below. 



the only diagrams of this type contributing to quantities 
of interest are those shown in Figs.|8]^c-f). We stress that 
all six diagrams in Fig. |8] contribute at the same order in 
WChPT. 

We begin by discussing the three-pion scattering di- 
agram (a). Only the w'^ vertex in £lo has the (tt'^)^ 
form needed to contribute. A straightforward calcula- 





(e) 



(f) 



FIG. 8. Tree-level diagrams in PQWChPT contributing to 
linear or quadratic dependence in Czw{t)- The flavor indices 
of the external fields are (implicitly) ordered as in Fig. [t] 
Notation for vertices as in Fig. [3] 



tion along the lines of those discussed in the main text 
leads to 



W-7 



1 



P L6(2Mo)3 



t + 



I 



3Afo 



(B5) 



The term linear in t arises, as usual, because the inter- 
action can occur at any time in the range — t. The 
volume suppression is now [compared to for 

the two-pion interaction, as in Eq. ([T])] because all three 
pious must be in contact. Note that at the order we are 
working in this appendix, Mq and M^r are interchange- 
able. 

Turning now to Fig. [sj^b), we find that only vertices 
having the pionic form (tt^) contribute. Such a form 
arises only from the mass, Wq and w'g terms in >Clo 



[eq. (15l]. The kinetic term does not contribute. After a 
straightforward but tedious calculation, we fincF] 



1 



t 

Mo 



2 
II 



L6(2Mo)4 

p-2Mot 



(B6) 



In the LCE regime, wj, ^ Mq//^, so the term linear in t 



is of the same order as that in the contribution (B5 1 



Recall that wL 



3/- 



and 



X 

7^' 



-2wL 
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The quadratic term in ( B6 1 arises because of the pres- 
ence of two vertices, much hke the quadratic term dis- 
cussed in the main text. The hnear term arises both 
from contributions in which the two vertices are close in 
time and integrated together from — t, and when one 
of the vertices is close to either or t. The latter origin 
suggests that one should also consider contributions in 
which one of the vertices is "absorbed" into the sources. 
Indeed, such diagrams, exemplified by those in Fig. [sjjc- 
f), do contribute at the same order. 



As is evident from the result ( B6 1 there are exponen- 
tially falling terms for small t. There are also contri- 
butions from excited states which show up in ChPT as 
contact terms. To avoid these, we consider henceforth 
only terms proportional to t and t'^. It turns out that 
there are no other LO diagrams leading to quadratic de- 
pendence, and the only other diagrams leading to linear 
dependence are those of Fig. [sj^c-f). Evaluating these 
diagrams, we find 



Higher order contributions will lead to corrections sup- 
pressed powers of M^/(47r/), a and 1/L. Assuming these 
corrections to be small, the coefficient of allows one to 
determine Wg , while that of t gives a combination of 
and Wg. Alternatively, w'g can be determined from other 
quantities such as the two pion correlators discussed in 
the main text. Either way, given w'^ one can use R^„{t) 
to determine Wy. 



A noteworthy feature of the result ( B8 1 is that the t 



term becomes comparable to the linear term at the rela- 
tively short time t ~ I/Mq. This differs from the corre- 
sponding results for the two-pion correlators [Eqs. ( |104[ ) 
and ( 105 1] for which the quadratic term becomes impor- 
tant at 



^(M„L)(/Lf»A_ 



(B9) 



6- ® 



Combining these results gives 



t + const. . 



(2Mo)4L6 



t 



(2Mo)3/2L6 



9 P 



(B7) 



(B8) 



[see Eq. ([5|]. The last inequality follows since one must 
have MqL ^ 1 and fL > 1 to avoid large finite- 
volume effects. The upshot of this discussion is that the 
quadratic term is much more important for the thrce-pion 
ratio than for the two-pion ratios. This is not true for 
cubic and higher order terms, which become important 
only at the longer times of Eq. ( B9 ) . 
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